Genomic attributes of airway commensal bacteria and mucosa.

Microbial communities at the airway mucosal barrier are conserved and highly ordered, in likelihood reflecting co-evolution with human host factors. Freed of selection to digest nutrients, the airway microbiome underpins cognate management of mucosal immunity and pathogen resistance. We show here the initial results of systematic culture and whole-genome sequencing of the thoracic airway bacteria, identifying 52 novel species amongst 126 organisms that constitute 75% of commensals typically present in heathy individuals. Clinically relevant genes encode antimicrobial synthesis, adhesion and biofilm formation, immune modulation, iron utilisation, nitrous oxide (NO) metabolism and sphingolipid signalling. Using whole-genome content we identify dysbiotic features that may influence asthma and chronic obstructive pulmonary disease. We match isolate gene content to transcripts and metabolites expressed late in airway epithelial differentiation, identifying pathways to sustain host interactions with microbiota. Our results provide a systematic basis for decrypting interactions between commensals, pathogens, and mucosa in lung diseases of global significance.

Crosstalk with lung fibroblasts shapes the growth and therapeutic response of mesothelioma cells.

Mesothelioma is an aggressive cancer of the mesothelial layer associated with an extensive fibrotic response. The latter is in large part mediated by cancer-associated fibroblasts which mediate tumour progression and poor prognosis. However, understanding of the crosstalk between cancer cells and fibroblasts in this disease is mostly lacking. Here, using co-cultures of patient-derived mesothelioma cell lines and lung fibroblasts, we demonstrate that fibroblast activation is a self-propagated process producing a fibrotic extracellular matrix (ECM) and triggering drug resistance in mesothelioma cells. Following characterisation of mesothelioma cells/fibroblasts signalling crosstalk, we identify several FDA-approved targeted therapies as far more potent than standard-of-care Cisplatin/Pemetrexed in ECM-embedded co-culture spheroid models. In particular, the SRC family kinase inhibitor, Saracatinib, extends overall survival well beyond standard-of-care in a mesothelioma genetically-engineered mouse model. In short, we lay the foundation for the rational design of novel therapeutic strategies targeting mesothelioma/fibroblast communication for the treatment of mesothelioma patients.

Molecular assessment of mycobacterial burden in the treatment of nontuberculous mycobacterial disease.

Introduction: Nontuberculous pulmonary disease causes significant morbidity and mortality. Efforts to tackle infections are hampered by the lack of reliable biomarkers for diagnosis, assessment and prognostication. The aim of this study was to develop molecular assays capable of identifying and quantifying multiple nontuberculous mycobacterial (NTM) species and to examine their utility in following individual patients' clinical courses. Methods: DNA was extracted from 410 sputum samples obtained longitudinally from a cohort of 38 patients who were commencing treatment for either Mycobacterium abscessus or Mycobacterium avium complex or who were patients with bronchiectasis who had never had positive cultures for mycobacteria. NTM quantification was performed with quantitative PCR assays developed in-house. Results: The molecular assays had high in vitro sensitivity and specificity for the detection and accurate quantification of NTM species. The assays successfully identified NTM DNA from human sputum samples (in vivo sensitivity: 0.86-0.87%; specificity: 0.62-0.95%; area under the curve: 0.74-0.92). A notable association between NTM copy number and treatment (Friedman ANOVA (df)=22.8 (3), p≤0.01 for M. abscessus treatment group) was also demonstrated. Conclusion: The quantitative PCR assays developed in this study provide affordable, real-time and rapid measurement of NTM burden, with significant implications for prompt management decisions.

M ycobacterium avium complex genomics and transmission in a London hospital.

BACKGROUND: Non-tuberculous mycobacteria (NTM) are environmental microorganisms and opportunistic pathogens in individuals with pre-existing lung conditions such as cystic fibrosis (CF) and non-CF bronchiectasis. While recent studies of Mycobacterium abscessus have identified transmission within single CF centres as well as nationally and globally, transmission of other NTM species is less well studied. METHODS: To investigate the potential for transmission of the Mycobacterium avium complex (MAC) we sequenced 996 isolates from 354 CF and non-CF patients at the Royal Brompton Hospital (London, UK; collected 2013-2016) and analysed them in a global context. Epidemiological links were identified from patient records. Previously published genomes were used to characterise global population structures. RESULTS: We identified putative transmission clusters in three MAC species, although few epidemiological links could be identified. For M. avium, lineages were largely limited to single countries, while for Mycobacterium chimaera, global transmission clusters previously associated with heater-cooler units (HCUs) were found. However, the immediate ancestor of the lineage causing the major HCU-associated outbreak was a lineage already circulating in patients. CONCLUSIONS: CF and non-CF patients shared transmission chains, although the lack of epidemiological links suggested that most transmission is indirect and may involve environmental intermediates or asymptomatic carriage in the wider population.

Evidence of immunometabolic dysregulation and airway dysbiosis in athletes susceptible to respiratory illness.

BACKGROUND: Respiratory tract infection (RTI) is a leading cause of training and in-competition time-loss in athlete health. The immune factors associated with RTI susceptibility remain unclear. In this study, we prospectively characterise host immune factors in elite athletes exhibiting RTI susceptibility. METHODS: Peripheral blood lymphocyte flow cytometry phenotyping and 16S rRNA microbial sequencing of oropharyngeal swabs was performed in a prospective elite athlete cohort study (n = 121). Mass cytometry, peripheral blood mononuclear cell (PBMC) stimulation and plasma metabolic profiling was performed in age-matched highly-susceptible (HS) athletes (≥4RTI in last 18 months) (n = 22) compared to non-susceptible (NS) (≤1RTI in last 18 months) (n = 23) athletes. Findings were compared to non-athletic healthy controls (HC) (n = 19). FINDINGS: Athletes (n = 121) had a reduced peripheral blood memory T regulatory cell compartment compared to HC (p = 0.02 (95%CI:0.1,1.0)) and reduced upper airway bacterial biomass compared to HC (p = 0.032, effect size r = 0.19). HS athletes (n = 22) had lower circulating memory T regulatory cells compared to NS (n = 23) athletes (p = 0.005 (95%CI:-1.5,-0.15)) and HC (p = 0.002 (95%CI:-1.9,-0.3) with PBMC microbial stimulation assays revealing a T-helper 2 skewed immune response compared to HC. Plasma metabolomic profiling showed differences in sphingolipid pathway metabolites (a class of lipids important in infection and inflammation regulation) in HS compared to NS athletes and HC, with sphingomyelin predictive of RTI infection susceptibility (p = 0.005). INTERPRETATION: Athletes susceptible to RTI have reduced circulating memory T regulatory cells, metabolic dysregulation of the sphingolipid pathway and evidence of upper airway bacterial dysbiosis. FUNDING: This study was funded by the English Institute of Sport (UK).

Sputum Proteomics in Nontuberculous Mycobacterial Lung Disease.

BACKGROUND: Nontuberculous mycobacterial (NTM) infections are difficult to diagnose and treat. Biomarkers to identify patients with active infection or at risk of disease progression would have clinical utility. Sputum is the most frequently used matrix for the diagnosis of NTM lung disease. RESEARCH QUESTION: Can sputum proteomics be used to identify NTM-associated inflammatory profiles in sputum? STUDY DESIGN AND METHODS: Patients with NTM lung disease and a matched cohort of patients with COPD, bronchiectasis (BE), and cystic fibrosis (CF) without NTM lung disease were enrolled from two hospitals in the United Kingdom. Liquid chromatography-tandem mass spectrometry was used to identify proteomic biomarkers associated with underlying diagnosis (COPD, BE, and CF), the presence of NTM lung disease defined according to American Thoracic Society/Infectious Diseases Society of America criteria, and severity of NTM. A subset of patients receiving guideline-concordant NTM treatment were studied to identify protein changes associated with treatment response. RESULTS: This study analyzed 95 sputum samples from 55 subjects (BE, n = 21; COPD, n = 19; CF, n = 15). Underlying disease and infection with Pseudomonas aeruginosa were the strongest drivers of sputum protein profiles. Comparing protein abundance in COPD, BE, and CF found that 12 proteins were upregulated in CF compared with COPD, including MPO, AZU1, CTSG, CAT, and RNASE3, with 21 proteins downregulated, including SCGB1A1, IGFBP2, SFTPB, GC, and CFD. Across CF, BE, and COPD, NTM infection (n = 15) was not associated with statistically significant differences in sputum protein profiles compared with those without NTM. Two proteins associated with iron chelation were significantly downregulated in severe NTM disease. NTM treatment was associated with heterogeneous changes in the sputum protein profile. Patients with NTM and a decrease in immune response proteins had a subjective symptomatic improvement. INTERPRETATION: Sputum proteomics identified candidate biomarkers of NTM severity and treatment response. However, underlying lung disease and typical bacterial pathogens such as P aeruginosa are also key determinants of the sputum proteomic profile.

Integrated genomics point to immune vulnerabilities in pleural mesothelioma.

Pleural mesothelioma is an aggressive malignancy with limited effective therapies. In order to identify therapeutic targets, we integrated SNP genotyping, sequencing and transcriptomics from tumours and low-passage patient-derived cells. Previously unrecognised deletions of SUFU locus (10q24.32), observed in 21% of 118 tumours, resulted in disordered expression of transcripts from Hedgehog pathways and the T-cell synapse including VISTA. Co-deletion of Interferon Type I genes and CDKN2A was present in half of tumours and was a predictor of poor survival. We also found previously unrecognised deletions in RB1 in 26% of cases and show sub-micromolar responses to downstream PLK1, CHEK1 and Aurora Kinase inhibitors in primary mesothelioma cells. Defects in Hippo pathways that included RASSF7 amplification and NF2 or LATS1/2 mutations were present in 50% of tumours and were accompanied by micromolar responses to the YAP1 inhibitor Verteporfin. Our results suggest new therapeutic avenues in mesothelioma and indicate targets and biomarkers for immunotherapy.

Airway microbial communities, smoking and asthma in a general population sample.

BACKGROUND: Normal airway microbial communities play a central role in respiratory health but are poorly characterized. Cigarette smoking is the dominant global environmental influence on lung function, and asthma has become the most prevalent chronic respiratory disease worldwide. Both conditions have major microbial components that are incompletely defined. METHODS: We investigated airway bacterial communities in a general population sample of 529 Australian adults. Posterior oropharyngeal swabs were analyzed by sequencing of the 16S rRNA gene. The microbiota were characterized according to their prevalence, abundance and network memberships. FINDINGS: The microbiota were similar across the general population, and were strongly organized into co-abundance networks. Smoking was associated with diversity loss, negative effects on abundant taxa, profound alterations to network structure and expansion of Streptococcus spp. By contrast, the asthmatic microbiota were selectively affected by an increase in Neisseria spp. and by reduced numbers of low abundance but prevalent organisms. INTERPRETATION: Our study shows that the healthy airway microbiota in this population were contained within a highly structured ecosystem, suggesting balanced relationships between the microbiome and human host factors. The marked abnormalities in smokers may contribute to chronic obstructive pulmonary disease (COPD) and lung cancer. The narrow spectrum of abnormalities in asthmatics encourages investigation of damaging and protective effects of specific bacteria. FUNDING: The study was funded by the Asmarley Trust and a Wellcome Joint Senior Investigator Award to WOCC and MFM (WT096964MA and WT097117MA). The Busselton Healthy Ageing Study is supported by the Government of Western Australia (Office of Science, Department of Health) the City of Busselton, and private donations.

Y disruption, autosomal hypomethylation and poor male lung cancer survival.

Lung cancer is the most frequent cause of cancer death worldwide. It affects more men than women, and men generally have worse survival outcomes. We compared gene co-expression networks in affected and unaffected lung tissue from 126 consecutive patients with Stage IA-IV lung cancer undergoing surgery with curative intent. We observed marked degradation of a sex-associated transcription network in tumour tissue. This disturbance, detected in 27.7% of male tumours in the discovery dataset and 27.3% of male tumours in a further 123-sample replication dataset, was coincident with partial losses of the Y chromosome and extensive autosomal DNA hypomethylation. Central to this network was the epigenetic modifier and regulator of sexually dimorphic gene expression, KDM5D. After accounting for prognostic and epidemiological covariates including stage and histology, male patients with tumour KDM5D deficiency showed a significantly increased risk of death (Hazard Ratio [HR] 3.80, 95% CI 1.40-10.3, P = 0.009). KDM5D deficiency was confirmed as a negative prognostic indicator in a further 1100 male lung tumours (HR 1.67, 95% CI 1.4-2.0, P = 1.2 × 10-10). Our findings identify tumour deficiency of KDM5D as a prognostic marker and credible mechanism underlying sex disparity in lung cancer.

A large-scale genome-wide association analysis of lung function in the Chinese population identifies novel loci and highlights shared genetic aetiology with obesity.

BACKGROUND: Lung function is a heritable complex phenotype with obesity being one of its important risk factors. However, knowledge of their shared genetic basis is limited. Most genome-wide association studies (GWASs) for lung function have been based on European populations, limiting the generalisability across populations. Large-scale lung function GWASs in other populations are lacking. METHODS: We included 100 285 subjects from the China Kadoorie Biobank (CKB). To identify novel loci for lung function, single-trait GWAS analyses were performed on forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC in the CKB. We then performed genome-wide cross-trait analysis between lung function and obesity traits (body mass index (BMI), BMI-adjusted waist-to-hip ratio and BMI-adjusted waist circumference) to investigate the shared genetic effects in the CKB. Finally, polygenic risk scores (PRSs) of lung function were developed in the CKB and their interaction with BMI's association on lung function were examined. We also conducted cross-trait analysis in parallel with the CKB using up to 457 756 subjects from the UK Biobank (UKB) for replication and investigation of ancestry-specific effects. RESULTS: We identified nine genome-wide significant novel loci for FEV1, six for FVC and three for FEV1/FVC in the CKB. FEV1 and FVC showed significant negative genetic correlation with obesity traits in both the CKB and UKB. Genetic loci shared between lung function and obesity traits highlighted important biological pathways, including cell proliferation, embryo, skeletal and tissue development, and regulation of gene expression. Mendelian randomisation analysis suggested significant negative causal effects of BMI on FEV1 and on FVC in both the CKB and UKB. Lung function PRSs significantly modified the effect of change in BMI on change in lung function during an average follow-up of 8 years. CONCLUSION: This large-scale GWAS of lung function identified novel loci and shared genetic aetiology between lung function and obesity. Change in BMI might affect change in lung function differently according to a subject's polygenic background. These findings may open new avenues for the development of molecular-targeted therapies for obesity and lung function improvement.

Bacterial Signatures of Paediatric Respiratory Disease: An Individual Participant Data Meta-Analysis.

Introduction: The airway microbiota has been linked to specific paediatric respiratory diseases, but studies are often small. It remains unclear whether particular bacteria are associated with a given disease, or if a more general, non-specific microbiota association with disease exists, as suggested for the gut. We investigated overarching patterns of bacterial association with acute and chronic paediatric respiratory disease in an individual participant data (IPD) meta-analysis of 16S rRNA gene sequences from published respiratory microbiota studies. Methods: We obtained raw microbiota data from public repositories or via communication with corresponding authors. Cross-sectional analyses of the paediatric (<18 years) microbiota in acute and chronic respiratory conditions, with >10 case subjects were included. Sequence data were processed using a uniform bioinformatics pipeline, removing a potentially substantial source of variation. Microbiota differences across diagnoses were assessed using alpha- and beta-diversity approaches, machine learning, and biomarker analyses. Results: We ultimately included 20 studies containing individual data from 2624 children. Disease was associated with lower bacterial diversity in nasal and lower airway samples and higher relative abundances of specific nasal taxa including Streptococcus and Haemophilus. Machine learning success in assigning samples to diagnostic groupings varied with anatomical site, with positive predictive value and sensitivity ranging from 43 to 100 and 8 to 99%, respectively. Conclusion: IPD meta-analysis of the respiratory microbiota across multiple diseases allowed identification of a non-specific disease association which cannot be recognised by studying a single disease. Whilst imperfect, machine learning offers promise as a potential additional tool to aid clinical diagnosis.

The fungal airway microbiome in cystic fibrosis and non-cystic fibrosis bronchiectasis.

BACKGROUND: The prevalence of fungal disease in cystic fibrosis (CF) and non-CF bronchiectasis is increasing and the clinical spectrum is widening. Poor sensitivity and a lack of standard diagnostic criteria renders interpretation of culture results challenging. In order to develop effective management strategies, a more accurate and comprehensive understanding of the airways fungal microbiome is required. The study aimed to use DNA sequences from sputum to assess the load and diversity of fungi in adults with CF and non-CF bronchiectasis. METHODS: Next generation sequencing of the ITS2 region was used to examine fungal community composition (n = 176) by disease and underlying clinical subgroups including allergic bronchopulmonary aspergillosis, chronic necrotizing pulmonary aspergillosis, non-tuberculous mycobacteria, and fungal bronchitis. Patients with no known active fungal disease were included as disease controls. RESULTS: ITS2 sequencing greatly increased the detection of fungi from sputum. In patients with CF fungal diversity was lower, while burden was higher than those with non-CF bronchiectasis. The most common operational taxonomic unit (OTU) in patients with CF was Candida parapsilosis (20.4%), whereas in non-CF bronchiectasis sputum Candida albicans (21.8%) was most common. CF patients with overt fungal bronchitis were dominated by Aspergillus spp., Exophiala spp., Candida parapsilosis or Scedosporium spp. CONCLUSION: This study provides a framework to more accurately characterize the extended spectrum of fungal airways diseases in adult suppurative lung diseases.

The undiagnosed disease burden associated with alpha-1 antitrypsin deficiency genotypes.

Alpha-1 antitrypsin deficiency (AATD), mainly due to the PI*ZZ genotype in SERPINA1, is one of the most common inherited diseases. Since it is associated with a high disease burden and partially prevented by smoking cessation, identification of PI*ZZ individuals through genotyping could improve health outcomes.We examined the frequency of the PI*ZZ genotype in individuals with and without diagnosed AATD from UK Biobank, and assessed the associations of the genotypes with clinical outcomes and mortality. A phenome-wide association study (PheWAS) was conducted to reveal disease associations with genotypes. A polygenic risk score (PRS) for forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) ratio was used to evaluate variable penetrance of PI*ZZ.Among 458 164 European-ancestry participants in UK Biobank, 140 had the PI*ZZ genotype and only nine (6.4%, 95% CI 3.4-11.7%) of them were diagnosed with AATD. Those with PI*ZZ had a substantially higher odds of COPD (OR 8.8, 95% CI 5.8-13.3), asthma (OR 2.0, 95% CI 1.4-3.0), bronchiectasis (OR 7.3, 95%CI 3.2-16.8), pneumonia (OR 2.7, 95% CI 1.5-4.9) and cirrhosis (OR 7.8, 95% CI 2.5-24.6) diagnoses and a higher hazard of mortality (2.4, 95% CI 1.2-4.6), compared to PI*MM (wildtype) (n=398 424). These associations were stronger among smokers. PheWAS demonstrated associations with increased odds of empyema, pneumothorax, cachexia, polycythaemia, aneurysm and pancreatitis. Polygenic risk score and PI*ZZ were independently associated with FEV1/FVC <0.7 (OR 1.4 per 1-sd change, 95% CI 1.4-1.5 and OR 4.5, 95% CI 3.0-6.9, respectively).The important underdiagnosis of AATD, whose outcomes are partially preventable through smoking cession, could be improved through genotype-guided diagnosis.

Presence of pleomorphic features but not growth patterns improves prognostic stratification of epithelioid malignant pleural mesothelioma by 2-tier nuclear grade.

AIMS: Nuclear grade has been recently validated as a powerful prognostic tool in epithelioid malignant pleural mesothelioma (E-MPM). In other studies histological parameters including pleomorphic features and growth patterns were also shown to exert prognostic impact. The primary aims of our study are (i) externally validate the prognostic role of pleomorphic features in E-MPM and (ii) investigate if evaluating growth pattern in addition to 2-tier nuclear grade improves prognostication. METHODS AND RESULTS: 614 consecutive cases of E-MPM from our institution over a period of 15 years were retrospectively reviewed, of which 51 showed pleomorphic features. E-MPM with pleomorphic features showed significantly worse overall survival compared to those without (5.4 versus 14.7 months). Tumours with predominantly micropapillary pattern showed the worst survival (6.2 months) followed by solid (10.5 months), microcystic (15.3 months), discohesive (16.1 months), trabecular (17.6 months) and tubulo-papillary (18.6 months). Sub-classification of growth patterns into high grade (solid, micropapillary) and low grade (all others) led to good separation of overall survival (10.5 versus 18.0 months) but did not predict survival independent of 2-tier nuclear grade. A composite score comprised of growth pattern and 2-tier nuclear grade did not improve prognostication compared with nuclear grade alone. Intra-tumoural heterogeneity in growth patterns is ubiquitous. CONCLUSIONS: Our findings support the incorporation of E-MPM with pleomorphic features in the epithelioid subtype as a highly aggressive variant distinct from 2-tier nuclear grade. E-MPM demonstrates extensive heterogeneity in growth pattern but its evaluation does not offer additional prognostic utility to 2-tier nuclear grade.

Utility of Nuclear Grading System in Epithelioid Malignant Pleural Mesothelioma in Biopsy-heavy Setting: An External Validation Study of 563 Cases.

Nuclear grading systems for epithelioid malignant pleural mesothelioma (MPM) have been proposed but it remains uncertain if they could be applied in a biopsy-heavy setting. Using the proposed system, we conducted an independent, external validation study using 563 consecutive cases of epithelioid MPM diagnosed at our institution between 2003 and 2017, of which 87% of patients underwent biopsies only. The median number of sites sampled was 1, with a median maximum tissue dimension of 17 mm (biopsy) and 150 mm (resection). The median overall survival (OS) was 14.7 months. The frequencies of grade I, II, and III tumors were 31% (132/563), 52% (292/563), and 17% (94/563). Grade I tumors were associated with the most favorable median OS (24.7 mo) followed by grades II (12.7 mo) and III (7.2 mo). The 2-tier nuclear grade separated tumors into low grade (19.3 mo) and high grade (8.9 mo). In multivariate analysis, 3-tier nuclear grade, 2-tier nuclear grade, and mitosis-necrosis score predicted OS independent of age, procedural type, solid-predominant growth pattern, necrosis, and atypical mitosis (all P<0.001 except 2-tier nuclear grade, P=0.001). In the scenario of a single- site biopsy with tissue dimension ≤10 mm, none but age (P=0.002) were independently predictive. Our data also suggested sampling 3 sites or a maximum tissue dimension of at least 20 mm from a single site is optimal for nuclear grade assessment. In conclusion our study confirmed the utility of nuclear grade in epithelioid MPM using a biopsy-heavy cohort provided the tissue sample met minimum dimensional criteria.

Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK Biobank.

BACKGROUND: Clinical and epidemiologic studies have shown that obesity is associated with asthma and that these associations differ by asthma subtype. Little is known about the shared genetic components between obesity and asthma. OBJECTIVE: We sought to identify shared genetic associations between obesity-related traits and asthma subtypes in adults. METHODS: A cross-trait genome-wide association study (GWAS) was performed using 457,822 subjects of European ancestry from the UK Biobank. Experimental evidence to support the role of genes significantly associated with both obesity-related traits and asthma through a GWAS was sought by using results from obese versus lean mouse RNA sequencing and RT-PCR experiments. RESULTS: We found a substantial positive genetic correlation between body mass index and later-onset asthma defined by asthma age of onset at 16 years or greater (Rg = 0.25, P = 9.56 × 10-22). Mendelian randomization analysis provided strong evidence in support of body mass index causally increasing asthma risk. Cross-trait meta-analysis identified 34 shared loci among 3 obesity-related traits and 2 asthma subtypes. GWAS functional analyses identified potential causal relationships between the shared loci and Genotype-Tissue Expression (GTEx) quantitative trait loci and shared immune- and cell differentiation-related pathways between obesity and asthma. Finally, RNA sequencing data from lungs of obese versus control mice found that 2 genes (acyl-coenzyme A oxidase-like [ACOXL] and myosin light chain 6 [MYL6]) from the cross-trait meta-analysis were differentially expressed, and these findings were validated by using RT-PCR in an independent set of mice. CONCLUSIONS: Our work identified shared genetic components between obesity-related traits and specific asthma subtypes, reinforcing the hypothesis that obesity causally increases the risk of asthma and identifying molecular pathways that might underlie both obesity and asthma.

Viral respiratory infections and the oropharyngeal bacterial microbiota in acutely wheezing children.

Acute viral wheeze in children is a major cause of hospitalisation and a major risk factor for the development of asthma. However, the role of the respiratory tract microbiome in the development of acute wheeze is unclear. To investigate whether severe wheezing episodes in children are associated with bacterial dysbiosis in the respiratory tract, oropharyngeal swabs were collected from 109 children with acute wheezing attending the only tertiary paediatric hospital in Perth, Australia. The bacterial community from these samples was explored using next generation sequencing and compared to samples from 75 non-wheezing controls. No significant difference in bacterial diversity was observed between samples from those with wheeze and healthy controls. Within the wheezing group, attendance at kindergarten or preschool was however, associated with increased bacterial diversity. Rhinovirus (RV) infection did not have a significant effect on bacterial community composition. A significant difference in bacterial richness was observed between children with RV-A and RV-C infection, however this is likely due to the differences in age group between the patient cohorts. The bacterial community within the oropharynx was found to be diverse and heterogeneous. Age and attendance at day care or kindergarten were important factors in driving bacterial diversity. However, wheeze and viral infection were not found to significantly relate to the bacterial community. Bacterial airway microbiome is highly variable in early life and its role in wheeze remains less clear than viral influences.